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§ 1. In this paper, without touching on the question of the existence of integrals 
of systems of simultaneous partial differential equations, I have given a method by 
which the problem of finding their complete primitives may be attacked. 

The cases discussed are two : that of a pair of equations of the first order in two 
dependent and two independent variables, and that of a single equation of the second 
order, with one dependent and two independent variables. 

I follow, as far as possible, the analogy of the method of Lagrange and Chakpit, 
and with this object introduce the conception of the " bidifferential " or differential 
element of the second order, which bears the same relation to a Jacobian taken with 
respect to two independent variables as a differential does to a differential coefficient. 

The solutions considered are, in general, complete primitives, that is, such as contain 
arbitrary constants in such number that the result of their elimination is the system 
of equations proposed for solution. The existence of such primitives is sufficiently 
established (see the papers of Frait von Kowalevsky and Professor Konigsberger, 
quoted hereafter) ; it will therefore be assumed, and the object of the investigation 
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will be to find conditions that must be satisfied by the equations of the solution and 
to put these conditions in a convenient form for solution by inspection. 

I should add that I am greatly indebted to the referees for their suggestions and 
for help in removing obscurities. 

To the list of authorities given by Dr. Forsyth (' Theory of Differential Equations/ 
Part L, pp. 299, 331), may be added the following : — 

Julius Konig, Math. Annalen, vol. 23, pp. 520, 521. 

Leo Konigsbergeb, Crelle, vol. 1.09, pp. 261-340. # Math. Annalen, vol. 41, 

pp. 260-28 5. f Math. Annalen, vol. 44, pp. 17-40. 
Ed. v. Weber. Mtinchen Ber., vol. 25, 423-442. 
J. M'Cowan. Edinb. Math. Soc. Proa, vol. 10, 63-70. 
Hamburger. Crelle, vol. 110, pp. 158-176. 
C Bourlet. Annales de l'ficole Normale (3), vol 8. 
Riquier. Comptes Rendus, vols. 114, 1.16, 119. Annales de l'Ecole Normale 

(3), vol. 10. 
Lloyd Tanner. Proc. Lond. Math. Soc, vols. 7-1 1. 
J. Brill. Quarterly Journal of Math., vol. 30, pp. 221-242. 

Several of the above papers are only known to me through abstracts. 

On BidifferentiaU. 

8 2. The idea of a " complete differential" plays an important part in the theory 
of differential equations. In this paper I shall try to show the importance of an 
extension of the same idea to differential elements of higher orders, such as enter 
into multiple integrals. 

An expression Xdx + Ydy is called a complete differential when X, Y are functions 
of the independent variables x 9 y 9 such that 

BY/dx = dX/dy. 

If this is the case, then, under certain restrictions, the value of \{Xdx + Ydy) depends 
only on the limiting values of the variables, and not on the intermediate ones by 
which these limits are connected, or, as generally expressed, on the path along which 
the integral is taken. 

This depends on the theorem that 



l(Xdx + Ydy) = f j (^£^)dxdy 

* For reasons stated below, I am not in agreement with the results given in the latter part of this 

paper. 

t In this paper it should be noticed that the equations (52) on p. 266 are not more general than (46). 
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when the single integral is taken round the boundary of the area over which the 
double integral is to extend. 

Further, X, Y are in this case the partial derivatives of a single function, 

§ 3. Let us consider the double integral 

| J(X dydz-{-Y dz dx + Z dx dy), 

where X, Y, Z are functions of the independent variables x, y, z. It is known 
that this, taken over a closed surface under certain restrictions, is equal to the triple 
integral 

1 1 fcdX/dx + dY/dy + dZ/dz) dx dy dz 

taken over the space enclosed by that surface. 

Hence, if dX./dx + dY/dy + dZ/dz = identically the double integral taken over 
a closed surface vanishes, and taken over two open surfaces with the same boundary 
has the same value ; that is to say, the value of the double integral depends on the 
values of x 9 y 9 z at the boundary only, and not, under certain restrictions, on the 
form of the surface enclosed by the boundary. 

By analogy we may call the element of the double integral a " complete double 
differential," or a " complete bidifferential " under these circumstances ; the condition 
that X dy dz + Y dz dx + Z dx dy may be a complete bidifferential is thus 

BX/dx + dY/dy + dZ/dz = 0. 

§ 4. A complete bidifferential may be expressed as a single term, such as du dv. 
For let u, v be two independent solutions of the equation 

dx dy dz 



so that u == a, v = b are integrals of the system 

cfec/X = dy/Y = dz/Z \ 

then x = *jj&£ Y = 0^4 z = e^, 

d(y, z) d(z, x) d(x, y 

6 being some multiplier, 

A 3X , 3Y , dZ d(0,u,v) 

and + + - 1 Aj-lJ 

ox dy dz d(x, y, z) 

Since the last vanishes identically 6 is a function of u 9 v only ; a function w of u, v 
may be found, such that dwjdu = 0, and thus 



d{y, z) y d(z, x) 9 l d(x, y) ' 
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Now in finding the value of the double integral taken over a part of any surface, 
it will be natural to suppose the co-ordinates of any point of such a surface to be 
functions of two parameters, say p, q 9 and to transform the integral into one taken 
with respect to these. The integral as transformed is 



if 



d(j q) d(p>q) B(p>q) m 



dp dq. 



and the known values in terms of p 9 q are to be substituted for x 9 y 9 z and their 
derivatives. 

The subject of integration is 

d(w 9 v) d(y, z) d(w> v) d(z, x) d(tv, v) d(x, y) 

d(y, z) d(p y q) d(z, x) d(p, q) d(x, y) d{ p y q)> 

dw dx , dw dy , dw dz dv dx , dv dy , 3^ dz 
die dp dy dp dz dp ' dx dp dy dp dz dp 



dw dx dtv dy dw d, 



dv dx _ dv dy , dv dz 



dx dq ] dy dq l dz dq ? dx dq r dy dq ~ r dz dq 
d(w t v) 



or 



d{p 7 q) ' 



The integral is therefore 



Cfd^v) , 
]}d(p } q) P q > 



dw dv 9 

6 



and if we take a single element we may write 

X dy dz + Y dz dx + Z dx dy = dw cfa;, 

dropping the parameters £>, q 9 since the values which x 9 ?/, 2; have in terms of them 
are immaterial. 

This equation is meaningless unless the expression in terms of parameters is under- 
stood. The same is true of ordinary differentials. If when u is a function of x 9 y 9 % 
we write 

7 dw 7 d-io 7 , du 7 

au = p— rfx + ~— a?/ + ?r~ oto, 

we mean that if a?, y, 2; are supposed to be any functions whatever of a single 
parameter p 9 then 

dw dn dx , du dy , 9^ cfe 
dp ^ dx dp dy dp dz dp * 

This equation being true quite independently of the expressions assumed for x, y 9 % in 
terms of p 9 we drop the denominator dp for convenience ; but in modern works on the 
Differential Calculus it is quite understood that a differential by itself is meaningless 
apart from this or some equivalent convention. 
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§ 5. The functions w, v are not uniquely determined. They may be replaced by 
W, V, where W, V are functions of. w 9 v 9 one of which, say W, is arbitrary, while V is 
only restricted by the condition 

9(W, V) 



o\w, v) 



1. 



The transformations of iv, v which are allowable will thus form a group. For a 
single integral the operations of the corresponding group consist in the addition of 
different constants, that is, in varying the constant of integration ; the theory of 
periodic functions is connected with discontinuous sub-groups of this. It is possible 
that an investigation of the discontinuous sub-groups of the group of transformations 
of two variables which leaves their bidifferential unchanged may lead to an extended 
theory of periodic functions of the two variables. 

§ 6. The finding of the functions tv, v may be considered as the indefinite integration 
of the bidifferential expression. It. is simplified by Jacobi's theory of the last 
multiplier, which is here a constant. 

KJ-LI.lv/t5 ii. " "■' r\. v • JL — — r\ , v 

6(y,z) d(z y x) 

we have X.dy — Ydx = ~~ div — tt dv \ 

cz cz 

and thus, on the supposition that v is constant, 

X dy — Y dx Y dz — Z dy Zdx — X. dz 



div 



dv dv dv 

dz dx xy 

_ 0? Z-v Y)dx + (yX-\Z )dy + (\Y - \x X)dz 

. dv , 9v , 9-y 
d£ dy cz 

Hence w may be found, if v is known, by integrating this last expression on the 
supposition that v is constant ; X, /x, v may have any values and the constant of 
integration is to be replaced by an arbitrary function of v. Thus, when one of the 
functions w, v is known, the other is found by ordinary integration. The only 
restriction on the one found first is the equation 

ex dy cz 

§ 7. Let us now suppose a greater number of independent variables. Let u be a 
function of x l9 x 2 ... x n . 
We have the relation 

cbw • — « djX\ j" ?\ aXp "t~" . . » "~v~ ^ Ci/Xft, 

CX i CX^ OOuii 



1 
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Here the differentials represent simultaneous infinitesimal increments, those of the 
ndependent variables being arbitrary. The equation may also be interpreted by 
supposing &*!, x % . . . x n to depend in any manner on a single parameter p, when the 

equation 

du ™ du dx r 
dp r= iOv- r dp 

holds whatever functions of the parameter we suppose x x . . . x n to be. 

To get the idea of a double differential we must suppose two sets of simultaneous 
infinitesimal increments ; denote them by d, 8. The bidifferential of x, y is then 
dx . 8y — 8x . dy* This vanishes if x, y are not functionally independent, just as dx 
vanishes if x is a constant. The analogy is very clearly shown if we say that dx 
vanishes when some function (\>{x) vanishes, dx dy vanishes when some function 
<j)(x, y) vanishes. 

If w, v are functions of n independent variables x i9 x % . . . x m we have 

CVU -— — 2Lf 77 ' CvXyy w ^*» 2L( ~— OX?* 

dv = 2 tt-; dx n 8v = 2 "^7 8a? r , and hence 

cfe . Si; — 8u . dv = 2 S ^7 ~- (cfav . Sa^ — efo* . S^ f ), 

or cfa^ tfc = 2 ^t^tt c^V d^ 

the summation being taken over all pairs of different suffixes r, s. Hence the 
expression for d^ dz; is formed by multiplying together 



with the conventions 



^ & , -j " dy , 

2^ r-v ClXf anCt 2*1 o> ~* OjXy 

j-> — J \JtX.',y if ~ \ \jJbf 



(ajJu (JjvJ — — """" U« Lt'.X', 
Cvt/y (a/Jl/ • — • U« 



We shall often use the notation c?(^, ?/) for efe t%. 

§ 8. For the purpose of double integration of such an expression as 2 X,, y d{x r> x s ), 

in which the coefficients X are functions of #j . . . sc # , it is natural to suppose x x . . . ^ 
expressed throughout the range of the integration in terms of two parameters, say 
p, q. The integral thus becomes 

* The dot is used here and throughout the paragraph to distinguish multiplication in the ordinary 
algebraic sense from multiplication according to the Grassmann conventions stated at the end of the 
paragraph. 
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X(t{J€ r} X s ji j j 
V9 7/ \ Ci/fJ (XQ • 

d(p, q) 1 1 



If X rs = or- 2 — ^ for all pairs of suffixes, the subject of integration in the last 

integral is d(u i v)/d(p, q), so that the integral becomes \\dudv. Its value will 
therefore only depend on the values of u, v 9 that is of x v x. 2 . . . x m at the boundary 
of the range of integration, and not on the form of the relations giving x Y , x 2 . . . 
in terms of p, q, which define the particular surface over which the integral is taken, 
In this case we may write 

r, s 

and call it a complete bidifferential. It is easily seen that the coefficients X satisfy 
the relations 

XrsXtj+XriXjs+XrjXsi^O, ....... (l) 

rVX" Pl"Y* r^Y 

3 1 ^y J rv *" """" v/« • • » . . . . » I^jl. 

t6| \JtX/g \JvUf 

for all combinations of suffixes, where it is understood that the term X rs d(x n x s ) may 
be also written X sr d(x s , x r ), so that 

^^ "Y" 

The conditions (2) are those which must be satisfied in order that the value of the 

double integral may depend only on the boundary. The difference of two values 

of the double integral, for which the same boundary is assumed, will be its value over 

a closed surface passing through the boundary curve, and this may be transformed 

into the triple integral 



\ 

si 



- i*OJ\. rs OJ\.- vr C-A..9J \ 7 7 j 

taken through the volume of any solid bounded by this closed surface. Hence this 
integral must vanish for any solid. By taking an infinitesimal solid, for every point 
of which all but x h x n x s are constant, we find the condition (2). 

The conditions (2) would be satisfied by an expression which was the sum of two 
or more complete bidifferentials, but (1) in general would not. 

§ 9. We next try to find whether these conditions are sufficient as well as 
necessary. Now all the coefficients X cannot vanish. Suppose that X 13 does not, 
then we have from (1) 

Y — ±2*j2*s ±1^±*Z („ « — q a v>\ 

and in virtue of these all the conditions (1) are satisfied. 
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Taking the values thus given for X m X,y, X^- we have 

aV !iY £*V 

m)= "5 h -5— + -5— 

7 d#$ cx s ox r 

^2s I @2\.] r . OJx.^ \ A lg /UxL^o. . C/xx.2* \ • -^*-lr f V-^-28 1 C^Vf: 

X 13 \ cte ?: cfc,, / X 12 \ o??i & f / X 13 \ 3^ 3^ 

+ -&-9/T I O &■ s \ . 02\-i^\ A ii / dxv. s 2 j 0-2v.2r\ -^2-/ /^Ms UA f | 

X 13 \ 0Xi Sx s j X 12 \ dx r dx s J X 13 \ dx r dx s 

Ap Cw£ ,A.-|q OXg ■"■"jo yJXf 

= f* (In) + |s (r2i) + |* (2si) + ^ (,li) + ^ (*2r) + ^ (l«r) 

A 12 ^12 ^13 A 12 -^13 A 12 



1 u 1 

\ ~/Y"Y_J_"\ r Y_l_"Y*'Y"\l _H. /"V* "V I V V 1 V "V \ 

+ f* (21s) + !* (21r) + |- (21i). 

A 12 A 12 A 12 

Thus the conditions (2) are not independent, but all follow from those in which at 
least one of the suffixes 1 , 2 enters. If they are satisfied then the equations 



n n 



can be satisfied by two integrals of the form w = a, 1; = b ; that is, these last 
equations will give x l9 x 2 as functions of the rest, such that 

For the conditions necessary and sufficient* for this are the vanishing of such 
expressions as 



J) Y V P} "Y" V Pi V PS -"V "V 3 V V^V 

1 v' b"" ^ 1 



,<?! 



cx s X 13 X 13 d^ X 13 X 13 cx 2 X 13 da? r X 13 X 13 r>x 1 X 13 X 13 3,2? 3 X 12 

in which 1 ? 2 may be interchanged and r, s are any two of the other suffixes. This 
expression may be written 

~ (rls) — —■ (12s) + vf" (12r) + ™ ^7 (X ]3 X„ + X rl X 2 « + X^ X ljf ), 

-A^ 3 -^-12 -^"12 xx j 2 Ot/q 

so that it vanishes and the conditions of integrability of the equations XX. lr dx r = 0, 
SX 2r daj r = are satisfied. If u = a, v = & are the integrals, then, since SX h /fav 
does not contain the differential of x l9 we must have 



^ 



For proof of this statement see Forsyth, 'Theory/ part I., pp. 43-51, 
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and 



In like manner 



Hence 



and 



2 X lr dx r = O i — civ 



r 



W 



dx x ^ dx x 



da 



jLi r — u 



-lr 



d(u, v ) 

(J [90*) Xy) 



Xo,. = 



-2r 



X 



rs 



(rst) 



6 A- 1 — - , the multiplier being the same. 

c($ 3 j x r ) 

9(0, u, v ) 

OiXf) 90$) v^t) 



Since this vanishes for all combinations of suffixes, 9 is a function of u, v, and if 
another function of them, w 9 is so chosen that 

dw/du = 6, 

we shall have 2 X^ d(x r> x s ) = 6 d(u, v) = d(iv> v). 

r, s 



Linear Differential Equations. 
10. If u — a is an integral of the linear partial differential equation 



*i 



>X, 



n + 1 



dx l 



-p X 



dx 



n+i 



•vy> GOG- 



% 



IXc 



4- MY 

i • • • i -*■*•# 



n+i 



dx, 



*&*n+u 



% 



where X l5 X 2 . . . X w+1 are functions of x u . . . aj w+1 , then ?6 satisfies the condition 



' T y — 



o, 



and the complete differential cfe is a linear combination* of the determinants 

O/X^) ax<£ 9 aXo » . . ax n9 ax^^^ 

the coefficients in the combination being usually functions of x l9 
If u = a is a common solution of the above equation and of 



« * . 



x. 



n+l* 






/ dx n + i 



f vX n j 



1 dx 8$ — -^w+ij 

then, in like manner, du is a linear combination of the determinants 

* This is generally expressed by saying that " a = a is an integral of the equations 



dx\ d%2 



tvtjLi 



n+1 » 



• » » 



XV " " * "V * 

1 "**-2 ■**-w -f- 1 

For the sake of the analogy with the work of §11, I prefer the phrase in the text, which expresses no 
more and no less than the one generally used. 
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Xxr -xr 

but in general, of course, It will not be possible to combine them so as to form a 
perfect differential. 

§ 11. An analogous process of integration may be given for two simultaneous 

equations 

£{Ay(p$ ; - — p/jti)} + tBipi + tdqi + E = 0' 

S{A'g(|)^ --##,■)} + SB'^ + SC'^ + E' = 0. 



(3), 



*».? 



in which the coefficients A, B, C, E, A', B', C, E' are functions of n independent 
variables, x l9 x% . » . #„, and two dependent y, 2, and 

% dz 

To fix the ideas, take n = 3 and let x ¥ x 6 stand for y % % respectively, A i4 for C;, 
A i5 for — B,-, A 45 for E, and make similar changes in the accented letters. Then, if 
u = a, v = b are two equations constituting a solution,* a, & being arbitrary con- 



stants, we must have 



"C^ 



2 A 

i, j = 1, 2 . . 5 






01 



C7( 1<V£, ^W/ 



i- 



V / 5 



*»J 



for 



1 l 3(y, ») 3fo, z) 



a, i? = &. 



and the values thus given for _p l9 g 1? j9 2 , j 2 , _p B , ^ 3 must satisfy the equations (3) identi- 
cally, since a, 6 are supposed arbitrary. The equations to be solved are thus reduced 
to others which are linear and homogeneous in the Jacobians, and which do not 
contain the dependent variables. 

The equations (4) give two of the Jacobians of u, v linearly in terms of the 
others ; if we substitute for these two in the identity 



d(u, v) = * ~^ A d(x h xj), 



vUii vV,;j 






we find that d(u 9 v) is a linear combination of the determinants of the matrix of tent 
columns. 



* This solution will not be a complete primitive unless a certain number of other arbitrary constants are 
involved as well as a, b, a supposition which is neither made nor excluded. 

It may be well to point out that the solution here assumed consists of two equations, and not of one 
equation involving an arbitrary function; in fact, any solution whatever necessarily consists of two 
equations, and one point of the present method is that these are to be sought together, not successively. 

t For n independent variables the number of columns in the matrix will be |(?i ■+• 1) (n + 2), the 
number of rows being still three. 
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(tMiX/'ij JUc) ) j CXl^ij Ct'o/j • 


> • Cvl «//^'j JUj J • 


, . CL\X> L ^ OC'^J 


A A 


A .. 


' • -"-45 


A' A' 

J*. 125 -**- 13? 


A'-- 


A' 

. XJL ^.g j 



There are thus eight bidifferential expressions, and the problem is to be solved by- 
finding such multiples of these as, when added together, will form a complete 
bidifferential. 

§ 12. As in the case of Lagrange's linear equation, this will generally, in practice, 
be done by inspection, and the method will be useful for finding solutions in finite 
terms — when such exist. But in any case, # whether the inspection is successful or 
not, there can be no doubt of the existence of suitable multipliers, in infinite number. 
For it is certain that the equations (3) have— possibly among other solutions — an 
infinity of solutions, each involving two arbitrary constants at least, and any one of 
these may be written u — a, v = b, where a, b are the two constants ; u, v are 
functions of the variables, but may, of course, be implicit functions of great com- 
plexity. The functions u> v must satisfy the conditions (4), and it immediately 
follows that d(u, v) must be a linear combination of the determinants of the matrix 
formed from (4) as above ; so that a corresponding system of multipliers must exist. 

If the solution is not in finite terms it is not likely to be found by inspection, and 
it is quite probable that the best way to find it would be by solving the original 
equations (3) in series. By whatever means the solution is found, the corresponding 
system of multipliers is thereby determined. 

If nine solutions of the form u = a, v = b have been found, the nine 
bidifferentials d(u l9 v x ) 9 diu^ v 2 ) . . . d(u 2i v 9 ) must satisfy identically a linear rela- 
tion, since they are all linear combinations of eight expressions only. 

We shall say that one of the nine pairs of functions is a " bifunction " of the other 
eight pairs. 

The following is, then, the definition of a bifunction. When the bidifferentials of 
any number of pairs of quantities are connected by an identical linear relation, with 
constant or variable coefficients, any one of these pairs is said to be a bifunction of 
the rest. 

The word bifunction is simply used as an abbreviation— at least for the present. I 
am not without hope that at a future time it may be found to have some connotation. 

* If one of the dependent variables with its derivatives is altogether absent from the equations (3), or 
if it can be made to disappear by a change of the other dependent variable, the equations (3) will in 
general have no solution. This case will then be excluded ; it is the only case in which the method of 
solution in series (as given, for instance, by Frau von Kowalevsky, ' Crelle/ vol. 80) cannot be used 
to prove that solutions actually exist. 

Another case that may fairly be excluded is that in which all the derivatives of one of the dependent 
variables do not occur or may be made to disappear by a change of the other. Such a system is equiva- 
lent to a single partial differential equation with one dependent variable, since the one whose derivatives 
are absent mav be eliminated. 

VOL. CXCV.— A. Y 
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It is, of course, evident that if u, v are functions of variables x v x 2 . . . then the 
pair u y v is a bif unction of all the pairs that can be formed from x l9 x 2 . . -. Other 
examples will be found later on in the paper. 

§ 13. Sometimes solutions exist for systems of partial differential equations in 
which the number of dependent variables is less than the number of equations. 

If, for instance, with the system just considered we take a third equation of the 
same form, the coefficients being distinguished by two dashes, there may be solutions 
common to the three equations. If u = a 9 v = h give such a solution, then it 
follows in like manner that d(u 9 v) is a linear combination of the determinants of the 
fo] lowing matrix : — 

Cv\X-\} Xy) » * . Cl[Xj 9 Xj) 

A A- 

A' A'» 






A " V 



Similarly for a greater number of equations. 

Application to other Differential Equations, 

§ 14. There are two classes of equations whose solution depends on that of a pair 
of such linear homogeneous equations as we have just been considering ; they are, 
firstly, systems of two equations in two dependent and tw r o independent variables, 
and, secondly, equations of the second order with one dependent variable and two 
independent. We shall consider them in order. 

Firstly, let y 9 % be the dependent variables and x l9 x z the independent ; sometimes 
we shall write x s for y and as 4 for z. Let p l9 p 2 be the partial derivatives of y and 
q l9 q% those of z 9 and let the equations be 

n\ x i> x z> V> % > .Pi' P& ft? ft?) = ®> 

J2\ x l> x 29 Vi Z > Pl> P& ft? ftO = o* 

A complete primitive will consist of two equations connecting x l9 x 29 y 9 % and 
involving four arbitrary constants. By differentiation these equations yield four 
more involving p l9 p Zi q l9 q 2 . As the two equations are supposed to be a complete 
primitive it must be possible to find expressions for the four arbitrary constants in 
terms of x v x 29 y 9 z, p l9 q l9 p 2 , q % ; the elimination of the four constants must give 

/i = 0,/ s =0. 

Let a l9 <%, a 3 , a 4 be the constants, and u l9 %, u^ u 4 the expressions for them m 
terms of x l9 x 29 y, z 9 p l9 q l9 p& q 2 . Suppose / 3 , f^f 59 f G to stand for u l9 u % , t/. s , % 
respectively. Then by differentiation we have for any value of the suffix i from 
1 to 6, 
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& + p & + 0r df} + ^ d?A + ^ ^ + ^ ^ + ^ -& = (r = 1 2) 

3^. •* r dy l -^ & 3p x <£r r 3p 3 &v 3^ 1 rfe r 3# 2 clx r ^ ' ^ ? 

the letter cl being used to denote differentiation with respect to x x or x 2 on the sup- 
position that the other is constant, while 3 indicates strictly partial differentiation. 

Since dp 2 /dx x = dpjdx^ dqjdx l = dqjdx 29 we find by eliminating the deriva- 
tives of p l9 q l9 p% 9 q% 9 that 

J(x l9 p l9 q l9 q 2 ) +p l J(y,p l , q }9 q 2 ) + q x J(z 9 p l9 q l9 &) + J{x 29 p 29 q ]9 q 2 ) 

+ .?%%> Pz> Vi> to) + ?2<% P*> 9i, &) = °> 
and J(x l9 q l9 p [9 p 2 ) + p x J(y 9 q l9 p l9 p 2 ) + q^(^ q l9 p l9 p 2 ) + J(%, £ 3 > Pi> P%) 

+ P2%> q»Pi> p%) + g&i(z, q2>Pi>.p*) = ° 

where J( ) denotes the Jacobian of any four of the functions f\> f%> f&f±>f&f§ with 
respect to the variables specified in the bracket. Of these equations there are thirty, 
but since they are given by the elimination of six quantities from twelve equations 
only six of the thirty can be independent. 

§ 15. One pair of these auxiliary equations will contain Jacobians of f l9 f 29 f 39 / 4 , 
and will in fact express the conditions that the equations 

dy = p l dx l + P>2 c l x z 

shall be integrable without restriction when p l9 p %9 q l9 q 2 have the values given by 
the equations /i = =f 29 / 3 = a l9 f Lh = <x 2 . 

Thus, if a pair of functions f. 69 J^ can be found satisfying these two auxiliary 
equations, the solution can be completed by solving a pair of simultaneous ordinary 
equations. (See Mayer's method, Forsyth, ' Theory of Differential Equations,' pp. 
59-62.) 

The two auxiliary equations that j^j f^ must satisfy are linear and homogeneous in 
their Jacobians, the coefficients of the Jacobians not involving the functions f^ 9 J\ ; 
the number of independent variables is apparently eight, but it may be taken as six, 
since two of the eight variables x l9 x 29 y 9 z 9 p [9 p 29 q l9 q 2 are given as functions of the 
other six by the relations f Y ^0,/ 3 =0, and may be supposed eliminated from f 39 f i9 
if that is desirable. 

The columns of the matrix formed as at § 11 are the rows of the following array :— 

d(x u y), 0, 0, 

d(x u z), 0, 0, 

d(x p y), 0, 0, 

(5) d(x p z), 0, 0, 

Y 2 
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c%, z), 0, S 

(10) d{z, p 1 ),q l ^q l , q 2 ], qi {p a , q,} + q,{p 2> q % }, 

d[%i,n), o, {qi,Pih 

d (y,p*),p-2{qi,q*}, pAqvPi}+pAq^Pi\, 

(15) dfo.ft). {q % ,p x }, {pvP%}> 

4*2. frMft.Pa}. 0, 

%l^l).Pl{fePll + Pa{<? 2 >.P 2 }> Pi{.Pi,P 2 J, 

^.ffi). qAq^Pi) + qAq^Pz\, qi{pi,p*}> 

c K x i>qz),{Pvqi}, o, 

d (y>&)>pAPi>vi} +Pz{p2>qi}>. pAp»pz}> 

d (^i> ffi), {»i»ff 2 } +JPi{y» ?a) + ftfc &}, -{*i,^ 2 ] -i>i{y,i> a } - ?i {«,#$}, 
(25) d{p lt q % ), - {x^qj - p x {y, q x \ -q&.qj, -{x z ,p 2 } -fofapz} -q z { z >P*}> 

d tPz, 9i),{»2. fe} + Pi{y,q*} + ftteft}, 1'^.Pi} + .Pi{y»2 J i} + ^iKp,}, 
#/i> &)>{»!. Pil +p 1 {?/,p i } + ^{^Pi} 4- {^,p 2 } +i?»{y» p 2 ) + qA z > &}, o> 

(5) 

Here {p l9 q l } 9 for instance, is written for d(f l ,f i )/d(p v q l ) 9 mi.d every fifth row is 

numbered. 

§ 16. In order, then, to solve the equations^ = 0, j^ = we have to form such a 

linear combination of the determinants of this array as will be a complete bidiffer- 

ential, say d(f B9 jQ, f 39 f± being such functions that the equations f Y = =f 29 f s = a 3 , 

.A = : c^can be solved for _p j? g 1? p 2f g 2 . The array contains twenty-eight rows, but 

thirteen of these are combinations of the other fifteen. For instance, multiply the first 

row by df } /dx 29 the second by dfjdy, the third by df } /dz 9 the seventh by dfjdpi> the 

eleventh by df } /dp 29 the fifteenth by dfjdq l9 the nineteenth by ^fi/dq^ and add ; the 

resulting row is 

cl(x l9 f } ) 9 0, 0, 

which vanishes. Other vanishing rows may be formed similarly by combining the 
rows of the array so as to have in the first column one of the following— 

<%i>/i)> d (x2,fi), <%,/i), <%/i), d (lhJi), d (thJi)> %i>/i)> d(&>fi)> 
d(x l9 fz) 9 d(x 2 J 2 ) 9 %,/ 2 ), d(z,f 2 ), d{p v f 2 ), d(p2,fz), d(q i9 f 2 ), d(q 2 ,f 2 ). 
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The coefficients in these combinations are partial derivatives of f Y or f^ thus, for 
instance, 

*(pl>A) = i x d (Pi> x i) + ^ d (Pi> x 2) + °iy d(pi> y) + {; d(pi, *) 

+ ^ d ('Pi> 9a) + lp 9 d (Pi> Pz) + ^l d (Pi> ( hX 

and so in other cases. 

The number of these combinations is sixteen, but it is to be lowered by three, 
since d(f v f l ) and d(f 2 ,f%) are identically zero and d(f l9 f 2 ) can be formed by com- 
bining the sixteen in two ways, so that three linear combinations of the sixteen 
bidifferentials vanish identically. 

Hence the array contains virtually only fifteen rows (28 — 13) and as there are 
three columns, we have thirteen bidifferential expressions to combine. Any pair of 
the four functions x l9 # 2 , y, z will satisfy the two auxiliary equations, as is clear either 
from the equations themselves or from an examination of the matrix ; of course these 
solutions of the auxiliary equations will not give a complete primitive. 

§ 17. If a complete primitive has been found it leads, as has been explained, to 
four equations 

and any pair of these must satisfy the auxiliary equations. Thus twelve pairs of 
functions satisfying these are known, namely 

Xi and Xj (i,j~ 1 ? 2, 3, 4) 

u t and Uj (i, j = 1, 2, 3, 4). 

These, however, are not all independent, but one pair is a bifunction of the other 
eleven. 

For if <£(#u x %> x & x ^ a i> c h> a 3> a i) = 

> 

y(^l5 ^25 "*%> ^4? ^1? ^2? ^3? a i) ™ ^ 

are the equations of the complete primitive, they must reduce to identities when 

i/ ]? u 2 , u 3 , u± are substituted for « 1? a 2 , a 3 , a 4i respectively. 
Hence 

<p\X\) &2> X Z-> X i? U l> %> ^3? ^4/ ' == ' ^ ] / . v 

XjjiX-t) Xc>) *£>?>) *^4j 1? / ^- / o? ^3? 4/ — ^ 

identically, and ^ 9<£> 7 _ ^ 3$ / 

i GXi i Ulli 

2* "^ OjXi — — *""■"" Jw ^ CLwiy 
■I O'Xj i Oil i 

so that 2 ^~— : d(x h Xj) = 2 x £%;, ?*,-) . . . . (8), 

and the bidifferentials of the twelve pairs of functions are connected by a linear 
relation. 



. . » . . » i 1 U I 
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§ 18. The method of Charpit for a single partial differential equation of the first 
order shows how all solutions may be deduced from one complete primitive, and it is 
a question of interest and importance whether there is any analogous method for 
simultaneous equations. Now it follows at once from the conditions for a complete 
bidifferential that a bifunction of the pairs that can be formed from m functions, say 
u l9 u % . . . . u m will be a pair of functions of w x . . . . u m . In the present case a 
bifunction of the six pairs that can be formed with u ly u %9 u 39 w 4 will be a pair of functions 
of these four, and the complete primitive to which it will lead will be the same as 
that given by u l9 u 2 . For when a solution of the auxiliary equations is known it leads 
directly to one and only one complete primitive by the integration of the equations # 

dy = Pl dx 1 + p 2 dx % 

. dz = q x dx x + q 9j dx% ; 

also the complete primitive to which the equations P 1 (u l9 u 2 , u 39 u±) = const., 
F 3 (u l9 u 9j9 u 39 u cl ) = const., will lead can be no other than is given by 

It must not, however, be forgotten that the system ¥ l = const., F 3 = const., 
f Y = 0, f 2 = may have a singular solution. If F 1? F 3 involve two other arbitrary 
constants this singular solution will involve four, and therefore in general be a com- 
plete primitive of the equations j^ = 0, f % = 0. Moreover, all new complete primitives 
are included among those thus given. 

For every solution implies six equations connecting x l9 x% 9 y 9 z 9 p l9 q l9 p 29 g 2 (two 
of these six are of course f x = 0, f % = 0), and, therefore, by . elimination of 
x l9 x 29 y 9 z 9 p l9 q u p% 9 q% 9 two equations or more connecting u l9 %i %9 u 39 u 49 which are 
known in terms of these eight quantities. If u l9 u 29 u 39 u^ are connected by four equa- 
tions they are constants, and the solution is therefore included in the old complete 
primitive. Let us, then, suppose that u^ a %9 t%, % are connected by two or by three 

equations, 

¥ a (u l9 %, u 3 , %) = (a = 1, 2 or 1, 2, 3). 

Now if pi, p%> g'u $%, are all expressed in terms of p 9 q 9 two of their number, and 
x l9 x %9 y 9 z 9 by means of the equations /[ = 0,/% = 0, the expressions 

dy — Pi^i "~" P%dXfy dz — q l dx i — ( l%dx % 
must both be expressible in the form 

k. x du Y + A 2 c£i% + A 3 rf% + A.^dti^ 

* Otherwise thus — -If in the auxiliary equations we suppose / 3 to have the known value u h they become 
a pair of linear equations for / 4 , which must be satisfied by u 2 , u^ u± ; now two linear equations in six 
independent variables can only have four functionally independent solutions, and one of these is known, 
namely, Wj . (In exceptional cases the two linear equations for u.>, w 3 , u± may be equivalent ; for instance, 
suppose /i = pi + qi, U\ = p 2 + &? fo having any form.) Hence, except in special cases, the particular 
complete primitive is defined when one of the functions u l9 u 2 , u Bl u. h or more generally a combination of 
them, F (ii v u^ 9 %, u±) is known, In the case supposed in the text two such combinations are known, 
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and since dp, dq are absent we must have in each case 



Thus the equations 



become 



4 a h- 

r = I <%> 



0, S A,. ^ = 0. 

d# 



4 

S 

r = 1 



d/?/ = p l dx l + p%dx<fr dz = q l dx l + g^cfotf 



2 



d/W-i) cLMoy ciUoy Ci/ii, 

di^ diic> dit% du it 

dp ' dp 9 dp 7 dp 

dit-^ du 2 du s du 4 

dq 9 dq p dq 9 dq 



4? 



0. 



These two equations, connecting du l9 du 29 du Z9 du^ taken with the system 



4 3F 

5 ^ du r = (a = 1, 2 or 1, 2, 3), 



= i du r 



show that if u ly u. 2 , ii s , ti 4i satisfy by themselves no other relations than P a = 
(a = 1, 2 or 1, 2, 3) we must have, as a consequence of the equations of the 
solution, 



i 3F. 



n 
cw 



.= i 9w,. dp 



4 3F a 3^., 

0, i -7T- -7T- = 0. 



3^. 3<7 



If, then, there are two equations 
the four equations 



F 



0, F 



2 



0, 



4 ^ ^ 

r = i 3i^ 3j9 



A 4 3F. du r 
r==1 era,, d^ 



(a = 1, 2) 



must reduce to two only. This will be the ordinary case, and we see that if the 
forms of F 1? F 2 , have been found by any means, the solution is completed without 
integration ; the process corresponds to Chabpit's method of deducing all complete 
primitives from one, but it differs in that the functions F l5 F 2 , are not arbitrary ; 
they must, in fact, be so chosen that the four equations last written shall reduce to 
two, and the conditions for this are clearly very complicated in general, though in 
particular cases available forms for F 1? F 3 may be seen on inspection. 
In the more uncommon case, when there are three equations 



the six equations 



x = 0, F 3 = 0, F 3 = ; 



^ 3W, 



i- = l Oil,, op 



4 d¥, "• 



°> ^S?^^" 1 ' 2 ' 3 )' 



must reduce to one only. 

These two cases are further discussed, from a somewhat different point of view, 
in SS 21—23. 

It should not be forgotten that the form in which the new complete primitive has 
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just appeared Is not that in which complete primitives were discussed in §14, since 
the equations are not here supposed to be solved for the arbitrary constants. 

§ 19. In addition to the six pairs (u i9 u/) of functions satisfying the auxiliary 
equations, we have also the six pairs (x i9 xj) ; of these twelve, eleve i are indepen- 
dent, the other being a bifunction of them. If we can find : bifunction of the 
eleven pairs which is not a bifunction of either set of six it will give a new complete 
primitive ; whether every, or indeed any, other primitive is thus given is a matter 
for further inquiry. 

Suppose Vi= hi (i = 1, 2, 3, 4) to be a new complete primitive, then it gives six 
more pairs of functions satisfying the auxiliary equations, and thus we have in all 
eighteen pairs. The bidifferentials of these must be connected by (18 — 13) five linear 
relations, one of which has been written (8) ; by means of the other four, an 
expression of either of the following forms— 

Ad(v l9 v z ) + Bd(v l9 v s ) + Gd(v l9 vj, 
Ad(v 29 t' 3 ) + "Bd(v z , v{) + Gd(v l9 1\ 2 ), 

can be found which will be equal to a linear combination of the twelve bidifferentials 
d (x i9 Xj) and d (u^ iij). It is natural to ask whether, conversely, any linear combina- 
tion of these twelve which can be written in one of the above forms will lead to a 
complete primitive ? In the first case this is not so, for if we take any function 
whatever, tj, of six independent variables, ij l9 g 3 . . . f G , we may choose the coefficients 



a l9 . . . a 6 , so that 

6 



2 ca d(r) 9 &) 



i-l 



shall be a linear combination of eleven* given bidifferentials ; the expression t a,; d^i 
may then be reduced to three terms, /? x <#£-, + /3 2 d^ + /3 3 <i£ 3 , so that for an 
arbitrary function (rj) a combination of the eleven given bidifferentials can be found 
of the ibvm fi l d(rj 9 £j) + /3 2 c%, £ 3 ) + & d{y) 9 Q 9 which is the same as Ad(v l9 v. 2 ) + 
Bd(v l9 v $ ) + Cd(v l9 t> 4 ). This argument does not apply to the second form 

Ad(v 29 v 3 ) + Bc% 3 , v { ) + Gd(v l9 v 2 ) 9 

and further investigation may show that any combination of the eleven that can be 
reduced to this formf will lead to a primitive. 

* Not of any lower number in general, since the most general bidifferential expression in this number 
of variables contains fifteen terms, while the expression just written vanishes identically if 

so that there are virtually only five coefficients, of which one must be left arbitrary. 

t The conditions necessary that a bidifferential expression may be reducible to this form include 
algebraic ones which are the same as for a complete bidifferential, since 

M(v 2i v B ) + Bd(v d9 vi) + Gd(vi,v 2 ) = -r-j A - dv ~ " ^ clVl } \-Adv* -Gdvi 
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§ 20. Before we can claim in any sense to have found the general solution of the 
auxiliary equations, we must be in possession of thirteen pairs of functions satisfying 
them ; we have only eleven when we know one complete primitive, and hence one 
more complete primitive, or even possibly two, must be found. An example (below, 
§29) will show that one more is not always enough. 

It is perhaps worth while to remark that any complete primitive defines the 
whole system of solutions, since it defines the differential equations. 

§ 21. The question of finding new solutions when a complete primitive is known 
may be attacked by the method of varying the parameters. Take the equations 
(6) or (7) of § 17. The problem is then to find such variable values for u u t% ? u $9 w 4 
as will satisfy the equations 



2 7r—dui = 0, 2 r^dui = (9). 



Since all variables are supposed functions of x l9 a%, we may make one of two 
suppositions with respect, to u l9 u% 9 u S9 u±\ either they are connected by three 
relations and are all functions of the same variable, say t 9 which is oi course a 
function of x l9 x %9 or they are only connected by two relations, so that two of them 
may be taken as functions of the other two. 

Suppose first that they are all functions of the one variable t. Then, generally, 
the four equations (7), (9) will define x l9 x %9 x< 69 x^ also as functions of t 9 and hence 
this supposition is not admissible unless it is possible to choose the functions of t in 
such a way # that the four equations (7), (9) will be only equivalent to three. The 

If these conditions are satisfied by an expression 

ij = 1,2 . . 6 

it can be put in the form 

2 Xi dxA ( 2 m dxi \ 

and then it must further be possible to express 

6 6 

2) Xi dxi and E /x^ dxi 

as linear combinations of three differentials, dv l9 dv 2 , dv B . The discussion of the conditions therefore 
belongs to the theory of the reduction of two such expressions, that is, of the extended Pfaff problem. 
* It seems obvious that this will not generally be possible ; but it may be well to give an example. 

Suppose the complete primitive to be 

y = axi 2 + bx-2 + c, 

z = cxi + ex-2 2 + hxiX-2 2 
so that the differential equations are 

V = iP i^i + j?2«2 + qi ~ i?2^2 2 , 

then the variations of the parameters a, b 9 c, e must satisfy the equations, 

xi 2 da + %%db + dc — 
Xidc + x-2 2 de + Xix 2 2 db = 0; 
VOL. CXCV.— A. Z 
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number of conditions, which will be of the nature of ordinary differential equations, 
thus imposed on the four parameters must not be greater than three ; for if they 
are subjected to four conditions they are made invariable ; it may be, however, 
less than three. For instance, a complete primitive of the equations Pi = p^ $i = % 
is given by 

y = a(x l + x 2 ) + b, z — c(x l + x%) + e ; 



he equations given by varying the parameters are 

• (x L + x 9 )da + db = 0, 
(x x + x 9 )dc + de = 0, 

which give the single differential equation connecting the parameters 

da de = db de. 

We may then assume arbitrary forms for two parameters in terms of a third, and 
find the fourth by integration. Say, for instance, 

b = <jf>(a), c = t/f(a), 
then e = \$\a)\\s\a)da, 

x i + % — — <£'( a ) ; 

thus we arrive at the known general solution 

V = x( x i + <>> * = w 0i + x ?)< 



whence, by elimination of Xi, 

(x^dh + dcfixtflb + cfc) + x^clehla = 0. 

This equation must fail to define x 2 , so that &, c, and a or e must be constant ; thence it follows that all 
four parameters must be constant. 

I lay stress on this, because it is not in agreement with the results of Professor Konigsberger (* Crelle," 
vol 109, p. 318), and appears in fact to show that his method there given is faulty. Professor Konigs- 
berger assumes (p. 313; I take m = 2) that the most general integral of the equations 

Mxi, %% y, z> pi> n> ft, #0 = 

has the form 

y = o>i(a;i, % 2 , h[h( x h ^)]» <fc#2<>4 5 %>]) 

z = (d 2 (^i, x 2 , <jf>i[^i(a;i, z 2 )], ^[^(^i, ^2)]), 

where <£i, c/> 2 denote arbitrary and ^ b ^ 2 definite functions. But suppose these equations solved for 
</>!, </> 2 in the form 

falfafa, Vi)] = X2^i, ^2, ?/, 3) 

and the arbitrary functions eliminated by differentiation. The differential equations thus formed are of 
the first degree in jh, lh, $u 2 2 > an( l are not D y an y means of the general form assumed. The differential 
equations in the examples given by Professor Konigsberger are, in fact, linear (see pp. 319, 328). The 
method appears to be founded on an interpretation of the last clause of § 2 (p. 290), which is not justified. 
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In the case of two equations of Clairaut's form 

V = Pi x i +^2% + <KPi> Pz> <Ii> %)> 
% = q lX[ + q 2 x 2 + xjj(p v pa, q l9 q 2 ), 

which will be more fully considered later, the number of differential relations among 

the parameters is two, so that one parameter may be taken as an arbitrary function 

of a second, and the other two found in terms of the second by solving two ordinary 

differential equations. 

If the primitive^ is 

y = a<% + b/3 + cy + <?S 

s = Aa + B£ + Cy + ES, 

where a, 6, c, e are the parameters, A, B, C, E known functions of a, b, c, e, and 
a, /?, y, 8 known functions of x i9 x 2 , then the variations of the parameters must 

satisfy the relations 

ada + fidb + ydc + Sde = 0, 

adA + (3dB + ydC + SeOE = 0, 

and thus, in general, if a, 6, c, e are all functions of one variable they are connected 

by three relations 

dA/da = dB/db = dC/dc = dE/de. 

The integral equivalent of these equations consists of three relations connecting 
a, 6, c, e with three arbitrary constants, and by eliminating a, 6, c, e we find a new 
solution of the original differential equations which is not a complete primitive, 
since it only contains three arbitrary constants. 

These examples show that the number of conditions to be fulfilled by the para- 
meters when all four are taken to be functions of one of them, may be one, two, or 
three ; this number is to be made up to three by assuming arbitrary relations (two, 
one, or none, as the case may be). 

§ 22. Usually the parameters will not be functions of one variable only, and we 
may suppose two of them, u s , u^ to be functions of the other two, u u u 2 . 

The partial differential coefficients 

du s du s du 4< du A 
dn Y ' dto 2 ' du l ' du 9j 

are then given by the equations (9), each of which is equivalent to two. The first, 

for instance, gives 

d<f> d<f> dn$ d<fi du 4i 

du x dtc s du x ^^// 4 da Y ' 

ty J dj> du 3 d(j> dv % __ 

3^2 du s du% du± du^ 

The derivatives are thus given in terms of ti v u. 29 w 3 , u^ x l9 x 2 , oj 3 , as 4 , and the last 

* It is unnecessary to give the differential equations, 

z 2 
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four may be eliminated by means of the relations (7) ; so that in the end we shall 
have two relations connecting u u %, %, %, and the derivatives ; the problem is of 
the same form as the original one, to solve two simultaneous partial differential 
equations in two dependent and two independent variables. 



Interchange of Variables and Parameters, 

§23. A curious thing may be noticed at this point. If in the equations <j> = 0, 
xjj = 0, we treat x Y , «%, x$, x Ai as arbitrary constants and eliminate them by differentia- 
tion, we are led to the same differential equations connecting u l9 a 2> u%, u^ as were 
just now given by the variation of parameters. Thus two equations in two sets of 
four quantities will give two pairs of simultaneous partial differential equations by 
taking each set of the quantities in turn as variables and the other as arbitrary con- 
stants. The auxiliary equations, if expressed in terms of the eight quantities, will, 
be the same in both cases ; this gives a meaning to the six solutions of the form 
(x h Xj) which we found the auxiliary equations to have, for any one of the six 
will lead to the primitive <£ = 0, \fj = of the second pair of differential equations, 
just as a solution (it,j, uj) leads to this primitive for the first pair ; any new solution 
of the auxiliary equations will in general lead to a new complete primitive for either 
pair, but an exception to this rule will arise when, for instance, the x differential 
equations have a complete primitive which gives three relations among u u u%, u^ t/ 4 . 

The array (5), transformed so that the variables are x ly x % , x 3 , x^ u ly %, w 3 , u^ 
connected by the equations <f> = 0, \fj = 0, will have six rows of the form 

six of the form d(u h uj), 0, 3 

and in the other sixteen there will be 

d(xi> Uj) in the first column, 



in the second the minor of 



3 3 <£ . 



dandify 



in the determinant 



d4 



Vf 



dxfiu^ ' dx } ' dx± 



3 3 <£ 



dx-idur, 

x & 

3 3 <£ 



dxjdti.^ 



dxcfiu-L ' ox 2 du 2 



^ 



d 2 cf> 

WXn\J'l'jr> 
li O 

o 3 </> 3 3 <£ 3 3 <^ 

dx$u l ' dx s du% ' dx g diin 

3 3 <£ ffij> 3 2 j> 

dx 4 du 2 ? dxjdu 3 



dxgdu h ' dx 2 

0(f) 



a 



X- 



'^0 



9 J 



dxodu 



3 3 c/> 



OXn 

o 



dx. 



dx^u x 



dcf> dcf> 

die 

3-^ 



dto 

d^fr 



d(f> d-\fr 
dxJSu^ ' dx % ' dx % 

3v 



(10) 



0, 



J. +j 



Ur 



du 



0, 



4 
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in the third the same expression with <f>, \jj interchanged. The array is thus practi- 
cally unchanged by interchanging the sets x and u 9 as should be the case. 
§ 24. This transformation may be accomplished by taking the equations 



i oui dx, i d'Uj ciXc 



y l * vw % c^ 2 



from which may be deduced 



d ( . 3$ dfwA c£_ / . 3$ ^ 



i' j dUiOUj cwx 2 dx l i au i ox l dx % $ oiiioy x cfe 3 $ 0^00 a^ 3 

_» ^ 3 3 $ *** rfwy ^ 3 3 (/> ^ ^ 3 3 <£ ^ ^ , ^ 3 3 </> dui 

i j dufi'iij d>x 1 dx% i duidx^dx-t ' i du- i dy J ^^dx l h i dii{dz^ dx^* 



ow p i9 q l9 p 2 , q$ are given by the relations 

tt~ 4- »i ^- + a-, V~ = 0, &c., 

cb x • - 1 l ay ll oz . 

and hence this equation, may be written 



2 



/fo? 3 dfa, y y z) 



% < 



dm \ dm 



CllX~t CtitA'Cfy <-f ) Zl 



( n ); 



in this <jf>, i/> may be interchanged so as to give another equation. 

Now, suppose 6 = a i9 x == % *° be two of the four equations connecting 
w l9 %, 1*3, ^ 4 with cc x , <% which yield a new complete primitive, and that y 9 z have 
been eliminated from 9 ^ by means of the equations <j> = 0, i// = 0, then the deriva- 

tives~T— \ ^p, &c, are given by the following relations : — 



9<£ diii 

Cflvi ChXi 



^ A. ' j~~ — ? 



. 3^ ^j 

^ _L _ — * — n 

v ??. ^ j ^L — 

and similarly for the derivatives with respect to # 2 . 

Substituting the values hence found for these derivatives in the equation (11), we 
have an equation linear in the Jacobians of the form 

Mai (i = 1, 2i ,■ = 1, *, s, 4), 
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the coefficient of the Jacobian written being the minor of %-• £7 in the determinant 



s t UWj 



(10). Hence the constituents in the second column of the transformed array are as 
stated, and those of the third are found in like manner. It is not, of course, neces- 
sary that these columns should be the same as would be found by actual substitution 
of the values of p l9 p 29 q i9 % in the columns of the original array ; a linear trans- 
formation is allowable, with constant or variable coefficients. 

The above process gives fifteen independent rows of the array ; the others are 
deduced from the consideration that y, z are known in terms of x v x 29 u l9 u 2i - u^ u^ 
from the equations <f> = 0, xjj = 0. 

Examples. 
§25. I. As a first example of the method of solution, take the equations 

a l — a 2> Pi = P& 

where a l9 f3 x denote known functions of x l} p l9 q x and a 2 , /3 2 known functions of 
x z> V& ft- 

In the array (5) multiply the seventh row by ~— — s , the fifteenth hy-^r^L the 

twenty-fourth by -g , *' \ - , and add. The result in the first column is d(a l9 /3 X ), in the 
second, by virtue of the particular forms of/j and f z , 

Kft) &i> ft} + &i> ft} { ( h>Pi} + (Pi, ft} {^i, ft} or 0, 
and in the third, 

{ x \>Pi) {p* ft} + {^i, ft} hh>P?} - {Pi> ft.} K> P2} or 0. 

Hence a 2 , ^ are two functions satisfying the auxiliary equations, and a solution is 
given by finding p j3 g 1? jp 2 , g 3 from the equations 

a x — a 3 = a, 

£1 = As = & 5 

and integrating. Two constants will be introduced by integration, so that the result 
is a complete primitive. 

§ 26. II. Take, secondly, the equations 

V = Pi»i + F (as 2 , i? 1? g 1? p 2 , &), 
2; = g 1 ir 1 -f Or (%, jp^ ft? P%> ft)- 
Here the twenty-fourth row is 

d (Pi> ft)> °> °> 
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so that jp 1? q l are two functions satisfying the auxiliary equations, and the integral is 
to be found by putting p Y = a, q } = &. Thus we have 

y — ax l = F (as 3 , a, 6, p 2 , g s ), 

— ^ = G (x 2 , a, 6, p 2 , <7 2 ), 

Or 7? rr: F (£ a, &, 17', £'), 

£ = G (£ a, 6, V ? 0> 
where g = x 29 yj ~ y -- ax l9 £ = z — &a? 1? 

These are ordinary differential equations, the solution of which will involve two 
new arbitrary constants and so constitute a complete primitive of the original 
equations. 

§ 27. III. The equations 

V = Pi x i +.P*Pz + <f>(Pi>Pz> Vi> &)> 

z = q x x x + q 2 x 2 + $(p l9 p 2 , ft, g a ), 

are of special interest, because more complete primitives than one can be found. The 
obvious solution is p } = a l9 p 2 = %, q x = &i, ^ = & 2? 

2; = 6^ + h 2 x 2 + i//(a x , a 2 , & 1? & 2 ). 

Suppose a l5 a 2 , & i? & 2 to be variable, but functions of one variable only— say a l9 then 
their variations must satisfy the relations 

cc L c/a x + ^2 ^ a 2 + dtj* = 0, 

^1 °^1 + ^2 f ^2 4" ^ := 0* 

These define x l9 x 2 , and, therefore, also y 9 z as functions of a l9 unless the determi- 
nants of the matrix 

da l9 da %9 d$ 

db l9 db 29 dxfj 

vanish ; it is necessary, then, that these determinants should vanish. Thus 
a l9 &J, a %9 b 2 are connected by two ordinary differential equations. We may assume 
any third relation connecting them at will ; suppose b x = F(a 1 ), F denoting an 
arbitrary function. Then by integration we may suppose a 2 , & 2 found in terms of a Y 
Also a Y is connected with x l9 x % by the relation 

X " + X * da[ + da, ~ °» 

so that a l9 b l9 a 2 , & 2 are all known in terms of x l9 »r 2 , and by substitution the values 
of y 9 % are found. 
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§ 28. The solution may be verified. We have taken b l = F(a 1 ), a known but 
arbitrary function of a x , and %, & 2 other functions of a ]? such that 



3^ d\jr 



d-^r da 2 



db 

fj_ 

da x 
dyfr db% 



da c 



F'K) : 3? 



dcL 



Then we have the further relations 



F'K) 



3</> d(j> 



L da 1 



"~r 



30 & 2 3(£ e£6 3 



i 



36 



j v J da 2 da ± db 2 da x 



*y + ** r* + ^ + d i n*d + f- r + 

" day d«j aOj v ' cta 3 <faj 



3<£ ^6 3 
36 3 da x 

d^r dh 



o, 



db 2 da x 



0,. 



which are of course not distinct. Also 



so that 



y = a l x 1 + a 3 % + </> (a 1? a 2 , 6 1? 6 2 ), 



JPi = c h + 






x i ~r % ^ i~ 5„ "+" 3„. 77 + ~ 7T + 



30 rf& £ 



& 2 30 

<&&! 3a x 3« 3 da x l db x da x ' db 2 da x 



= a 



h 



and in like manner p 2 = a 



2' 



Again 2; = & x crr^ -j- 6 2 x % + *K#i> %? &i> &<&), and 



3: 



&i + 



da x 

CvdLs-1 



db 



Jbi 



da x 



1 + X 



clb 2 d-\Jr d^jr da% dyjr db x dyjr db 2 
h da x oa x da% da x Sb x da x d\ da x 



K 



and similarly q % = 6 2 . 



Hence the original -differential equations are actually satisfied. If the arbitrary 
relation assumed — which may if convenient involve more than two of the parameters- 
contains two arbitrary constants, the new solution will generally be a complete 
primitive, since two more constants are introduced by integration.* 



•x- 



The ordinary equations to be integrated may have a singular solution with one arbitrary constant, 
or with none : if the arbitrary function has been chosen so as to involve three or four arbitrary constants, 
the whole number being thus raised to four, the solution so given may quite well be a complete primitive, 
and, in general, will be so. 

t The above investigation in a modified form shows how to find integrals of a system of three 
equations 

fa (U, % pi, $% 01, £2) - 0, "N 



> 



5 » a 



a * 



• (12) 



fa (%v,l>hP2, qi, 22) - 0, 
fa (u,v,z>i,p2>qi, 02)-= o,J 

where u = jh %i + 1H %2 ~ y, v = q 1 % x + q 2 x 2 — z« 
One solution is to take ti, v, 2 J h 2 7 2, qi, ^2 as constants connected by the three relations (12) ; if they are 



not constants we have 



du = %i dpi + ^2 '^23 
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§ 29. Let us now consider the new solutions of the auxiliary equations, given by 
the new complete primitive. The old solutions are the six pairs of the form x h Xj and 
the six of the form u h Uj 9 where u x = p i9 u 2 = p 2 , u 3 == q l9 u_ % = q 2 . The bi- 
differentials of these twelve satisfy the relation 



d(y, z) — p } d(x. v z) -- p%d(x 2 , z) 



Xl + * ~ ■ a * 






mJh 



»i + 



4- 



&o + 



5 ? **a-> 


<K 


d\]r 


3t^ 


Spi 


<¥z 


3</> 


d^> 


¥? 


3?2 




■ 8? 3 


3<£ 


30 


3ft' 


3fe 


3ilr 


4. ^ 

3? 3 



rf (Pi.2>a) + 



ai + 3K' 



3?i 
3^ 



3ft' ai + 3ft 



d iPi, 9i) 



rf (Pi> %) + 



Xi ~j~ 



3-^ 



3£ 

3^r 



tjC 1- .l -m i 

dpi " l dqi 



d(p s , q Y ) 



d (Pv to) + 



30 

3ft' 



dcf> 



d(q lt q 2 ) 



«i + g£> % + ^ 



In the auxiliary equations we may take x u x 29 p l9 p 2 , q }9 q 2 as independent variables, 
since y, z are given explicitly in terms of these six. 

From (12) follow three more relations connecting the six differentials du, dv, dp h dp 2 , dqi, dq 2 , so that 

their ratios are determinate, and therefore u, v,pi, qi,p 2 , #2 can only be functions of one variable. The two 

equations last written will then, generally, give Xi, x 2 in terms of this variable, which may not be. Hence 

we must have 

du = Xdv, dpi — Xdqi, dp 2 -■ Xdq 2 , 

and since df\ - 0, df 2 = 0, df B = 0, and du, dv, dp h dq h dp 2 , dq 2 do not vanish, A must satisfy the 
equation : 



du * dv dpi dqi ' dp 2 ' dq 2 

X d ll + d k , \ d ll + d ll , A,_# + ^2 



\ d k + c Mi, x c M^ + ^ 

r/w ^ ? r/px dqi 



\ d fA + d Il 

dp 2 dq 2 



= 0. 



If X satisfies this equation the differential relations du — Xdv, dpi 
since u, v, p h q h p 2 , q 2 are connected by the equations 



Xdq h dp-2 =• A^ 2 reduce to two only 



h - o,/ 2 = o,/ 3 



0. 



By integrating these two we find two more relations involving two arbitrary constants. Hence we may 
suppose v, pi, p 2 , q h q 2 expressed in terms of u, and -find a solution by eliminating u from the following : — 

/ u « pix x + p 2 x 2 - y, 
v = qi xi + q 2 x 2 - z, 
1 = Xj dpi/du 4- x 2 dp 2 jdu. 
VOL. .CXCV. — A. 2 A 
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Then 
d(x l9 y) = pA x i> •%) + ( x i + fyJ d ( x i>Pi) + ( x z + fyj d ( x i> Pa) + 9^ d ( x u 9i) + 

g^T d \ x l> $2/9 

and similar expressions may be found for d(x 29 y), d{x v z) 9 d(x 29 z) in terms of the 
bidifferentials of the pairs of independent variables. 

Let c }9 /* 2 , c 3 , c^ be the constants of integration in a new complete primitive 
found by the method of §§ 27-8. Let X be the common value of the ratios 
dp ]9 /dq l9 dp 2 /dq 29 d<f>/d\fs. Then, after integrating the equations dpjdq x = 
dpjdq 2 = dcj)/dx(j (= X) by help of an assumed relation connecting, say, p 1? g l5 p l9 g 3 , 
e 1? c 3 we have four relations among 

Pv P& 2lJ fe ^? ^1? C 2? #8 6 4? 

and we may therefore suppose p l9 p %9 q l9 q 2 expressed in terms of X, c l9 c 29 c 3 , c 4 , 
unless X is a constant, and therefore itself a function of c 1? e 29 c 39 c 4 . Then 

cZpj — X cfaft , c?p 3 — X c?g. 3 , d<j) — X dxjj 

will be linear combinations of dc l9 dc 29 dc %9 dc^ and so will some such expression as 

a dp x + fi dX 9 

where a vanishes if X is one of the constants or a function of them. Conversely, 
dc l9 dc 29 do ?)9 dc^ will be linear combinations of 

dp } — - X dq l9 dp 2 — X dq 29 d<j> — X d\f* 9 ot dp l + fi c£X, 

and the bidifferentials of e l? c 2 , c 3 , c 4 in pairs will be linear combinations of the six 
following expressions : — • 

d (Pi> Pa) - Xd (?i> Pa) - XcZ (jPi» 2a) + X * d (?i> 2a)> 

^ Cpi> <£) — ^ (g^ <£) — X< ^ ( Pi> VO + x ^ (<Zi> */>)> 

d (Pa> <£) — X ^ (#a> <£) — Xc? (Pa> VO + ^ (fe */>)> 
fid (p l9 X) — aX^ (q [9 p Y ) — /3X(i (g l9 X), 

ad(p 29 p Y ) + £d(p 2 , X) — aXd(q 29 pj — fiXd(q 29 X), 

ac£ (<£, p x ) + fid (<jf>, X) — aX<i (i|f, ^j) — /3Xc£ (i/r, X). 

These are combinations of the bidifferentials of p 1? p 29 q l9 q 2> in pairs, with the 

expressions 

d (p l9 X) ~ Xc? ($ l9 X), 

d {p 29 X) - Xd (q 29 X), 

<i (<£, X) — X<i (t//, X). 
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Now X is a definite function of x x , a; a) p x , p 2 , q x , q%, given by eliminating the 
differentials from the equations 

dp x = \dq x , dp % = hdq z , d(j> = kd\}/, 
By means of the first two, the third becomes 



dp x dg x 



oft oft 



m«7i + n- + ^ - \ 2 -* - x^ 

f ' 3Pg fyg ^ft <*ft 



dg 3 = 0, 



and the fourth 



v 4. X ^t 4- ^ 



><% + {^3 + X g^ + fy W = ° 



The result of elimination is therefore 



1 \ ^2 9 & 



X 3 ^ 



- X 



Bl|r V 



dp 2 dq, 

9Pa <W 



x^ + ^Vx— + — Wo 

Oft <W\ &I Oft 



This shows the form of X as a function of x x , x%, p x , p%, q x , q 2 , not involving c x , c 3 , 
c 3 , c 4 . Now this choice of X makes it possible to choose coefficients A, B, C, E, F, G, 
such that 



x x dp x + x % dp z + d<f> = A (c^jOj — \dq x ) + B (c£p 3 — X<% 2 ) + C (d(f> — \d\jj), 
x i fy.1 + x 2 ^a + ^ = E (dp x — \dq x ) + F (<ip 3 — Xc?^ 3 ) + G (d<f> — Xdi/>). 



Thus 



A{rf(2? ls X) - Xd(q x , X)} +B {d(pz, X) - Xc% 2 , X)} 

+ G{d(cj), X) — Xe£(^, X)} = x l d(p 1 , X) + x z d(p 2 , X) + cZ(<£, X) 

= multiples of bidifferentials of p x , p 2> q x> q 2 



+ 



+ 



3\ 
dx x 



x x d (p x , x x ) + x % d (p z> x x ) + d (<j), x x ) 
x x d (p x , x % ) + x % d (p %> x, 2 ) + d ((j>, x. 2 ) 



+ multiples of bidifferentials of p } , p%, q x , q % . 

2 A 2 
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In like manner 

E {d( Pl , X) - U{q x , X)} + F {d{ lh , X) - \d(q. 2 , X)} 

+ G {d (^X) - Xd (<|/, X)} = ^ {d (z, x k ) - g. z d (x s , Xj)} + d ~{d {z,x„) - q Y d (x it %)} 

+ multiples of bidifferentials of_p l5 p 3? ft, ft. 

Hence the three expressions 

d(p l9 X) - \d(q l9 X), 

d((f) 9 X) — Xc£(i/f, X), 

are all reduced to the same, save for a factor, by adding or subtracting multiples of 
the bidifferentials of x l9 x %9 x Z9 x± and of u l9 %, u Z9 u^ ; the same is therefore true of 
the bidifferentials of c l9 c 2 , c 3 , c 4 . Hence all the new complete primitives found by 
the method of §§ 27-8 only add one to the eleven known " bifunctionally " indepen- 
dent pairs of functions satisfying the auxiliary equations ; one more pair, leading to 
a fresh complete primitive, is yet to be found. 

§ 30. These results may be used to construct examples of bifunctions. For 
instance, the equations 

y = pi x i + p% x % + 2i> 
2 = STl^l + % x % + Pz> 

lead to the following case among others :— 
In the equations 

d Il = ^ _ ^i _ x 
rfft dy 3 djpj 

put ft = X + a, tfft = c£X, and integrate. 

Thus 2p a = X s + 6, 3^ = X 3 + c, 4p x = X 4 + e, 

and the arbitrary constants a, &, c, 6 in the new solution are respectively equal to 

ft — X, 2p 2 — -X 2 , 3ft — X 3 , 4p x — X 4 , where \ 2 x x + # a + X = 0. 

Now from § 29 it follows that d (c, e) can be expressed in terms of c? (a, 6), the 
bidifferentials of x l9 x %9 y 9 z and those of p l9 p%, q l9 ft. 
For convenience, let us write 

u 9 v 9 to, x 9 y 9 z for x l9 X, p l9 p 29 ft, ft respectively ; then 
for x % we must put — * v(l + uv) 9 
for y „ „ wu ~ xv(l + w) + y, 

for z „ „ yu — zv(l + uv) + x 9 
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so that the eight original variables .connected by two equations are now expressed in 
terms of six. 

Thus d (3y — v 3 9 iw — v*) can be expressed in terms of d (z — v 9 2x — v*) 9 the six 
bifferentials of iv 9 x 9 y 9 z and those of 

%i 9 — v(l + uv), wa — xv(l + uv) + jy> yu — ^(1 + uv) + X, 

that is, of 

u 9 v 9 wu — aw(l + w) + 2/ 5 .yw — ■ zv(\ + w) + #. 

There is no difficulty in finding the relation. It is 

u*d(3y — v* 9 iw — v*) — 6v 2 (l + uvfd(z — v 9 2x — v 2 ) 

— 12uU(y 9 w) + 12v 2 (l + wv) s d(«, a;) 

+ 12?; 2 {(1 + uv) (y — zv*) — ?^(w — aw 2 )}d(i;, ^) 

— ■ 12^ 3 (1 + uv)d{v 9 yu — zv — uzv 9, + '#) 

+ 12wv s d(v, tra — aw — waw s + 2/) == 0. 

Here then we have an identical linear relation connecting the bidifferentials of 
seven pairs of functions of six variables. Any one of the seven pairs is accordingly 
by definition a bifunction of the other six. 



Second Application. 

§ 31. Take now a differential equation of the second order, 

ffay,z,p 9 q 9 r 9 s 9 t) = 0, 

where p 9 q are the first and r, s 9 t the second partial derivatives of z with respect 

to x 9 y» 

A complete primitive will consist of a single equation in x 9 y 9 z involving five 
arbitrary constants, say a l9 a 29 a 3 , <% a B . If we form the first and second derivatives 
of this equation we shall have, in all, six equations from which a l9 a% 9 a 39 a 4 , a 5 can 
be found in terms of x 9 y 9 z 9 p 9 q 9 r 9 s 9 t 9 and the original differential equation will be 
the result of eliminating a l9 a %9 a 3 , a 4 , a 6 . Let u i9 u 29 u 39 u^ 9 u 5 represent the 
expressions found for a l9 a 29 a 39 a 4 , a 5 respectively, in terms of x 9 y 9 z 9 p 9 q 9 r 9 s 9 t. 

Then from the equations 

j = 0, u } = a l9 % = a 29 

by differentiating, we can form six equations which will involve the third derivatives 
of z ; by eliminating these we deduce the following two differential equations to be 
satisfied by u l9 u % : — 
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J(x, r 9 1) + pJ(z, r, t) + rJ(p, r, t) + sJ(q, r, t) + J(y, s, t) + gJ(3» 5, 
+ sJ(p, s 9 1) + tJ(q, s 9 t) = ? and 

Here J ( ) denotes the Jacobian of /, w 1? w a with respect to the variables specified. 

These equations express the conditions which are necessary and sufficient in order 

that 

dz = pdx + qdy, 

dp = rdfo? + $dy, 

dq = sdte + £<% 

may be integrable without restriction, when r, s, £ are given in terms of x 9 y i z, p 9 q 9 
by the equations 

J- = 0, u x = a 1? ^ == % 5 

the conditions must of course be satisfied by any three of the six functions %, %, 
%> ^4> ^55 ./• We thus have forty equations, of which only eight can be algebraically 
independent. 

§ 32. The conditions to be satisfied by u l9 u % are linear and homogeneous in their 
Jacobians with respect to the eight variables x 9 y, z 9 p 9 q 9 r 9 s, t ; of these, one is 
given in terms of the rest by the equation f— 0, and may, if convenient, be sup- 
posed not to occur in u v u % : hence the auxiliary equations in this case have seven 
independent variables and the dependent variables do not occur explicitly : to 
find a solution we are therefore to form a complete bidifierential, which shall be a 
linear combination of the determinants of the following array :— 



d(r, s), 







— ■ A. - 


- j?Z — rP — sQ 


d(r,t), 


X -j- pZ -\- T P + 


S\£ 9 


V 


- gZ — sP — £Q 


d(s, t), 


Y + qZ + sP + ^Q, 







d(p,r), 


rT, 






- rS - sT 


(5) d(p,s), 


O JL. % 






rR 


d(p, t), 


— rR — • sS 






sR 


d(q, r), 


O JL. m 






*"™™" OlO """"*" v JL 


d(q, s), 


tT, 






sR 


d(q, t), 


"""*" SXaj —■ " 0k) 9 






tR 


(10) d(z, r), 


P T, 




- 


- pS — gT 


CtyZ 9 SJ 9 


q% 






^>R 


CviZ 9 0\ 9 


— pR — qB, 






gR 


d(x 9 r) 9 


T, 








Cf'iXy Si} 


0, 






R 
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(15) d(x 9 t) 9 - R 

%,r), 0, ~T 

d(y,s) 9 T, 

<$(?/, t) 9 ■— S, R 

d(a;,_p), 0, 

(20) d(y 9 p) 9 0, 

d(z 9 p) 9 0, 

d(x 9 q) 9 0, 

%,?), 0, 

d(z 9 q) 9 0, 

(25) d(p 9 q) 9 9 

<#(#,£), 0, 

%,*), . o, . 

d(x 9 y) 9 0, 

X, P . . . are written for dj/dx 9 df/dp . . . 

Of these twenty-eight rows, only twenty-one are independent. For instance, 
multiply the 1st, 2nd, 4th, 7th, 10th, 13th, 16th by — - S, — T, P, Q, Z, X, Y respec- 
tively and add ; the resulting row is 

d(f 9 r) 9 0, 0, 

which vanishes since f= by hypothesis. 

Suppose d (u l9 Uc>) to be the complete bidifferential formed from the determinants 
of the array, then to complete the solution we have to find r, s 9 t from the 

equations 

j •=. 0, u-± = a l9 u% = <% 2 ) 

and integrate the equations 

dz = pdx + qdy ? dp = rdx + sdy 9 dq = sdx + £<%. 

It will amount to the same thing if we treat u x as known in the auxiliary 
equations. They must be satisfied if u& u^ %i h are substituted in turn for u 2 . Now 
two homogeneous linear partial differential equations in seven independent variables 
can at most have five common solutions, and here one of these, u l9 is known ; the 
other four may be taken as ?i 2 , u Z9 n^ 9 u 5 . 

§33. Any two of the five functions x 9 y 9 z 9 p 9 q will satisfy the auxiliary equations, 
but as we have to solve for r, s 9 t, these solutions will not serve our purpose. They 
are ten in number, and ten more will be given by taking in pairs the expressions u l9 
u 29 u 39 u A9 u b given by any complete primitive. These twenty are not all bifunc- 
tionally independent, for since there are three relations^ among the ten expressions 

X > V> %9 P> 9> U l> ^2? ^3> ^4> ^5> 



-* 



Compare § 34, p. 184. 
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three linear relations can be formed connecting the twenty bidifferentials ; one is 
formed from each pair of equations as at § 17 (8), Hence seventeen bifunctionally 
independent solutions of the auxiliary equations are known when we have one com- 
plete primitive. The full number is nineteen ( -'- — 2 j, and in order to know all 

we must have one, or possibly two (see § 41, p. 190), more complete primitives. 

§ 34, New solutions found by varying the parameters may be divided into two 
classes, according as the parameters are or are not all functions of one variable ; 
solutions of the former class only occur in exceptional cases, and the principles of § 21 
apply to them with slight modification. 

Let the three equations connecting 

^> y> %> p.i !?> %* %j u & ^4» %> 

be fa(x, y, z 9 p, q, u l9 u %9 u& u^ u b ) = (i = 1, 2, 3) ; 

(the forms fa, fa 9 fa are not unrestricted, but must be such that the following rela- 
tions hold identically 

or we may take fa as not involving p 9 q and fa as p dfa/dz + hfafix 

fa as q dfa/dz +■ dfa/dy), 

then the variations of the parameters must satisfy the three equations 

2 ~^du r = (i = 1, 2, 3), 

in order that the same relations may subsist among x, y 9 z, p 9 q 9 r 9 s 9 t and the para- 
meters, as held when the parameters were constant. 

If the parameters are functions of one variable, their forms must be so chosen that 
the three equations last written reduce to one only, otherwise we shall have five 
relations connecting x 9 y 9 z 9 p 9 q with this single variable. 

§ 35. If the parameters are not functions of one variable, only the equations 

2 3 Z - du r = 

are equivalent to six, and determine the partial derivatives of u 3 , v %9 u B with respect 
to u l9 xi % in terms of the five parameters and x 9 y 9 z 9 p 9 q. By help of the relations 
fa = we may suppose x 9 y 9 z 9 p, q eliminated and thus arrive at a system of four 
partial differential equations connecting u l9 u %9 u^ u i9 %, 
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The original system may also be taken to cpnsist of four equations connecting five 
variables x 9 y, z, p, q 9 namely : 

dzjdx = p, dz/dy = q, dp/dy = dq/dx 
f(x, y 9 2, p, q, dp/dx, dp/dy, dq/dy) = 0, 

and so the method of variation of parameters does not lead to any simplification of 
the problem in general. 

§ 36. The interchange of variables and parameters is again possible ; it is, perhaps, 
made clearer by taking three equations of perfectly general form, 

(pi {X l9 X% 9 X* 3 , X^, &* 5 , U± 9 U% U§ 9 u^, u 5 ) = (i = 1, 2, o), 

connecting two sets, each of five quantities. 

Whichever set we' suppose constant and eliminated by differentiation, we are led 
to a system of four partial differential equations connecting the quantities of the 
other set, two of the five being taken as independent variables. A new solution of 
either of these systems of differential equations will in general yield a new solution of 
the other. 

Suppose, for instance, that we have a new solution of the u equations ; this gives 
u z , u^ u 5 , say, in terms of u l9 u 2 . Then the six equations included in 

S* 5^ du r — (i = 1, 2, 3) 

r=1 Cu r v 7 

give two relations among x l9 ... u l9 u %9 since the four differential equations, 
which are consequences of these six, are supposed satisfied ; by the help of these 
two, u i9 u^ 9 may be eliminated from the three relations <^» 1 = 0, <£ 3 ~ ®> $3 = 0> an( i 
thus three relations are given connecting x l9 x 2 , x 39 x±, x 5 ; these three will constitute 
a solution of the x system of differential equations. 

§ 37. In this more general case there will not seemingly, as a rule, be any more 
solutions for either system of differential equations. For the derivatives, say, of 
x Z9 x A9 x 5 with respect to x l9 x 2 are given in terms of these five variables and two 
others, say u lf -u^. The forms we may assign to u l9 u % are then restricted by three 
differential equations derived from the three conditions 

® Xr = fer . ( T = 3^ 4 5\ 

and thus, generally speaking, no forms of u l9 u 2 will be suitable. In some cases the 
conditions are not inconsistent, and we may form an array by the method of § 11 
such that if d(9 9 y) is a combination of its determinants, then 6 = a, x = b, <£ x = 0. 

VOL. CXCV. — A. 2 B 
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<jE> 3 = 0, 3 = will give suitable values for u u % This array will have four columns 
and forty-five rows, ten such as 

d (x i9 Xj) t 0, 3 3 
ten such as d (ui, uj) y 0, 0, 0, 

and twenty-five of the following type. In the first column there is d{x iy uj), in the 
(r + l)th the minor of &<j> r /dXiduj, in the determinant 



for 



a% 


d 3 -<j> r 


3^,. 


3 2 0>- 


3 3 0, 


S*i 


3*2 303 


dx Y du^ 


dx % du^ 
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Ml 

3c,j 

.V 




^1 

3z 5 








3^3 


M 

dx q 


3^3 


9 3 

3^4, 
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Ma 

dx y 


Ma 

dx 2 


303 

3^ 3 


303 

3^4 


303 











r 



X j Lj ? O. 



This interchange of variables and parameters may take place whenever their 
numbers are equal, the differential equations being of the first degree. 



Examples. 

§ 38. I. As an example of the method of solution take the equation a = /J, where 
a is a function of r, s i p — sy, cc and /3 a function of s, t 9 q — sx 9 y. 

In the array (§ 32) multiply the first row by 3a/3r ? the fifth by da/dp, the fourteenth 
by 3a/3^, the seventeenth by — s da/dp, and add ; the resulting row is 

d(a, s), 0, 0, 
Hence we take a = /3 = a, 5 = 6, 

2 = 6x^/ + X + Y, 

X being a function of x only and Y a function of y only. Then a = a is a relation 
connecting as, dX/dx, d 2 X/dx 2 , and /3 = a is a relation connecting y, dY/dy, d?Yjdy^, 
and by solving these for X, Y respectively we shall have the complete primitive. 
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§ 39. II. As a second example take the equation 

F (r, s 9 t 9 p — sy 9 q — ty 9 z — qy + J % 2 , x) = 0. 

Here the third row of the array is 

d(s 9 t) 9 0, 0, 

so that the functions s 9 t satisfy the auxiliary equations. Put, then, s = a 9 t = b ; 

thus 

q = ax -\- by -\- c 

% — axy + I by 2 + cy + X, 

the last term being a function of # only. The differential equation thus becomes 

F (d' 2> X./dx 2 9 a 9 b 9 dX./dx 9 ax + c 9 X, x) = 0, 

an ordinary equation of the second order giving X in terms of x and two more 
arbitrary constants ; hence the finding of a complete primitive is reduced to the solu- 
tion of the equation last written. 

§ 40. III. If the equation is of the particular form F(r, s 9 t 9 p — rx — sy 9 q — sx — ty 9 
z — - px —* qy-^-^rx^-^sxy-j-^ty 2 ) = 0, the first three rows of the array are 



d (r, s) 








d(r, t) 








d(s, t) 





0. 



Hence any two of the three functions r 9 s 9 t will satisfy the auxiliary equations, 
and a complete primitive is given by putting 

r = a 9 s = h 9 t = b. 

Hence p = ax-\-hy-\-g 9 q = hx~\-by-\-f 

z — c + gx +fy + J (ra 2 + 2to?/ + % 2 )> 

where a, &, c, jf, </, /^ are constants satisfying the relation 

F (a, h 9 b 9 g 9 f 9 c) = 0. 

This is a case in which other solutions are readily given by supposing the para- 
meters variable and functions of one variable only, say a. The variations must 
satisfy the conditions 

xhla + 2xydh + iy 2 <i& + 2xdg + 2ydf--\- 2dc = 0, 

xda -\- y dh -\- dg ■= 9 x dh -\- y db -{- df = 0, 
2 B 2 
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whence follows xdg-\-ydf-\-2dc =: 0, a simpler relation that may be taken Instead of 
the first of the three. 

These equations will define x 9 y in terms of the single variable a 9 unless all the 
first minors of 

da dh dg vanish. 

-if 77 7 p 

ah ah af 
dg df 2dc 

We thus have three ordinary differential equations connecting a 9 b 9 c 9 f 9 g, h ; 
they are connected also by the relation F (a, h, b 9 g, f 9 c) = 0, and the fifth relation 
among them may be chosen arbitrarily, so that we may put h = <f> (a), an arbitrary 
function. 

Then we have 

db/da = {ft (a)} 2 , df/da = ft (a) dg/da, 



2dc/da = (dg/da) 



2 



F (a, <£(a) 5 &, gr, /, c) = 

as the equations determining b 9 g 9 f 9 c in terms of a. These are to be integrated, 
and then a is to be eliminated from the equations 

x + y dh/da + dg/da = 0, 
z — c + ^as+j^ + i ( ra 3 -j- 2^7/ + % 3 ). 

The result of elimination will be a solution of the differential equation. Three 
constants are introduced by integration, and thus, If the function <j> involves two 
constants, the new solution will generally be a complete primitive. 

§ 41. The new complete primitive gives new solutions of the auxiliary equations 
which we shall now examine. Let a ]? a 3 , a 3J a 4 , a 5 be the new set of parameters. 
Then a, h, g 9 &, c, / are connected with these parameters by five equations, one of 
which is the original equation F = 0. These five relations are such, that if 



then 



dh = Xda 9 dg = /xcfe , 
c?& = \ % da 9 2dc = /A£a, cZf = \/*c£a ; 



of these five, the first two define X, jjl in terms of a, a 1? a.,, a 3 , a 4 , a 5 , and the others 
must then follow from the five equations that give /&, #, &, c, / in terms of a and 
the same new constants. Thus, In general, we may suppose a, h 9 g 9 &, c 9 f 9 /x, 
expressed in terms of X, a 1? a 2 , a 8 , a 4 , a 3 and the expressions will be such that 



dh — Xc/a, <% — /xc?a, rf& — X s c£a, 2c?c — jt*, s c£a, c^/ 1 — Xfida 
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involve only the differentials of a 1? a 3 , a 3 , a 4 , a 5 . One of these five is expressible in 
terms of the other four, since 



3F 3F 3F 3F 3F 

da da + dh dh + Zb dh + d g d 9 + W df+ ^ de= °' 
while one of the relations connecting A, p, a, b, . . . is 

3F , JF , 3F , x9 3F , . 3F , , n 3F 

3« + X dh + ^37 + X 3T + ^3/ + ^37 = °" 

Some expression such as vdk — pda will also involve the differentials of a T , a 2 , a 3 , 
a 4 , a 3 only. Hence the differentials of a ]9 a a , a 3 , a 4 , a 5 will be linear combinations of 
vd\ — pcfa, <iA — • Xrfa, (% — ^cda, db — \ 2 da, df — \u,da, 2dc — /x 2 cfa, of which the 
last five satisfy a linear relation. 

Thus the bidifferentials of a l9 a 2 , a 3> a 4 , a 5 in pairs will be linear combinations of 
the bidifferentials of a, &, c, /, $, h (only five of the six need be used) in pairs, and 
of the expressions 

d(h 9 X) — \d(a, X), d(g, X) — [id(a 9 X), d(&, X) — X 2 (f (a. X), 
d(/, X) — \jjid(a 9 X), 2d(c, X) — fL 2 d(a, X), 

of which last five, only four are independent. 

Now X, /x are connected not only by the equation 

3F , x 3F , 3F , X9 3F , x 3F , , 9 3F 
37+ X 3T+^37 +X 3T+^37+^ 3 37= ' 

but also by the equation 

x + \y + p = 0, 

so that they are definite functions of a?, y-, a, 6, /, g 9 h. 

Again p = aas + % + g 9 

d(p 9 x) — M(y, a?) = xd(a, x) + yc?(&, x) + c?(g^, a;), 
d(p, y) - ad(aj, 2/) = xd(a, y) + yc?(&, #) + d(g 9 y). 

Thus y{d(h 9 X) — Xc£(a, X)} + {d(g 9 X) — /x<i(a, X)} 

= xd(a 9 X) + yd(/&, X) + d(g 9 X) 

- d~ i d (p> x ) ~ hd (y> x )l + %y l d (p> y) - ad ( x > v)\ 

+ multiples of bidifferentials of a, b, c, f, g, h. 
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In like manner 

y{d(b 9 X) — X^d(a 9 X)} + {d(f 9 ^) — hpd(a, X)} -f* x{d(h, X) — \d(a, X)j 
= xd(h 9 X) + yd(b 9 X) + d(J 9 X) 

+ multiples of bidifferentials of a, &, c, /, gr, A e 
Lastly, 

2xy{d(h, X) — \d(a, X)} + y*{d(b 9 X) — X 3 d(a, X)} + 2x{d(g 9 X) — /j,d(a, X)} 

+ %y{d(f, X) — Xfid(a, X)} + (2cZ(c, X) — jj?d{a 9 X)} 

= a? 2 c?(a, X) + 2xyd(h, X) + y*d(b 9 X) + 2xd(g 9 X) + 2yd(f 9 X) + 2<i(e, X) 

= 2{d(z, x) — (lix + % + f)d(y 9 x)}d\/dx 
+ 2{d(z 9 y) — {ax + % + g)d(w, y)} d\/dy, 

+ multiples of bidifferentials of a, 6, c, /, (7, A. 

Hence, in all, nine combinations of the ten bidifferentials of a l3 a %9 a 3 , a 4 , a 5 can 
be expressed in terms of the bidifferentials of x 9 y 9 % > p 9 q and of a, 5, c, j^ (/, A ; that 
is, in terms of the bidifferentials of the seventeen known independent pairs of functions 
satisfying the auxiliary equations : thus the new complete primitive adds only one to 
the number of these known bifunctionally independent pairs, and one more must be 
added in order to give the full number. 

This theory enables us again to construct examples of bifunctions of a number of 
known pairs which may reach eighteen. 

§ 42. The foregoing investigation may be modified so as to give singular solutions 
of a pair of differential equations of the form in question, say 



Fi(n s, t 9 p 9 q 9 z) = 

F s (r, s 9 t 9 p, q 9 z) = 



where p = p — rx — sy, 

q = q — sx — ty 9 



/*9 *-* ■' ■■ — ■ zy 



i(p + p)x - ±(q + q)y. 



A complete primitive would be given by supposing r, s 9 t 9 p 9 q 9 z constants con- 
nected by the above equations. Another solution would be given by solving the 
total differential equations found by supposing the relations 
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x dr + y ds + dp == 0, 
as c?s + y dt + rfg = 0, 

x dp + y cZg + 2cfe = 0, 

to reduce to the same relation linear in x and y. That is, we must solve the system 

ds = Xc/r, <i£ = X 2 cZr, dp = jac£r, 
where X, /a are given in terms of p, g, z, r, #, t by the relations 

dr' + x aT + x "aT + ^aj" + V ^ + *^ ar - o, 
aT + x aT + x ar + **$ + x ^ + ^V = ° 5 

and ^, 2 in terms ol p, r, 5, £ by the relations F x = 0, F 2 =0. 

The complete primitive of these ordinary equations will involve three arbitrary 
constants, and there may be singular solutions with a lower number ; none of these 
will therefore constitute a complete primitive of the partial differential system 

F x = 0, F 3 = 0. 



